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Abstract 

The dynamics of Rydberg states of atomic hydrogen illuminated by reso- 
nant elliptically polarized microwaves is investigated both semiclassically and 
quantum mechanically in a simplified two-dimensional model of an atom. 
Semiclassical predictions for quasienergies of the system are found to be in 
a very good agreement with exact quantum data enabling a classification of 
possible types of motion and their dynamics with the change of the elliptic- 
ity of the microwaves. Particular attention is paid to the dynamics of the 
nonspreading wave packet states which are found to exist for an arbitrary 
microwave polarization. 
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I. INTRODUCTION 



Consider a hydrogen atom, initially in a high Rydberg state with the principal quantum 
number n Q illuminated by the microwave field of a frequency uj close to the frequency of 
the unperturbed Kepler motion ujk = l/ n o- Quantum mechanically speaking, a resonant 
periodic field couples strongly several n states due to almost resonant transitions n — > n' = 
n ± 1 for n close to n . Since the driving field is periodic, applying the Floquet theorem 
U one may find eigenstates of the atom-field system (the so called Floquet or dressed |2| 
states). The eigenenergies are then referred to as quasienergies of the system and are defined 
modulo hu. The Floquet states, periodic in time, may be viewed as linear combinations 
of unperturbed system eigenstates. While the construction of Floquet states is possible 
both for nonresonant and resonant driving, in the latter case they may have quite unusual 
properties, especially in a semiclassical limit. 

Classically a resonance between the driving frequency and the frequency of the unper- 
turbed motion leads to a strong perturbation of the system and a creation of a stable island 
in the phase space centered on a periodic orbit of the frequency uj. The motion in the island 
is locked to the microwave frequency due to the nonlinear resonance ||. Semiclassically one 
then expects that the corresponding Floquet time-periodic state will follow a classical tra- 
jectory (in the vicinity of the periodic orbit) i.e. form a wavepacket which will not disperse 
in time. 

States localized in the resonance island for such a periodic perturbation have been first 
considered more than 20 years ago [|J and details of their semiclassical construction for some 
one-dimensional (ID) model systems have been analyzed fj|||. The wavepacket character of 
the time evolution of individual Floquet states has been realized only quite recently for 
hydrogen atom driven by linearly polarized microwaves. Independently, it has been shown 
that Gaussian wave packets may propagate almost without dispersion along circular periodic 
orbits in hydrogen atoms driven by circularly polarized fields || . The fact that the harmonic 
approximation implied in M and resulting the Gaussian wave packet form is not a necessary 



condition for non-spreading properties have been discussed in PJTOll where it was shown 
that the harmonic expansion provides a good approximation for exact Floquet states of the 
system (which by definition, being time periodic, do not spread on a long time scale). 

The fact that one may construct in a nonlinear system wavepackets that do not spread 
induced a flurry of activity in the field. Some [p7T| p~6[] concentrated on modifying the potential 
so as to make the harmonic approximation as good as possible aiming at the construction 
of Gaussian nonspreading wavepackets. Claiming that the anharmonic corrections are big 
for a hydrogen atom driven by circularly polarized fields those authors added a magnetic 
field perpendicular to the microwave polarization plane []IT|--|TE|]. This helps to minimize the 
unharmonicities in this plane leaving, however, unaffected the terms along the magnetic 
field axis. Thus the Gaussian wave packets remain even then an approximation to the real 
dynamics and must disperse (although very slowly in time). Another approach aimed at 
optimizing the coordinate system to the symmetry of the problem [fL7| . 

As discussed by us elsewhere P, |IT||i8| -f2~T[| much more fruitful is another approach, already 
outlined above. Namely we define the nonspreading wavepacket as a single Floquet state 
(for which the Gaussian packet may be merely an approximation) . Then a localization of the 
wavepacket in the vicinity of a stable fixed point is assured by the correspondence principle 
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provided the size of the surrounding island in the phase space is comparable to %. One 
may construct an approximate resonance Hamiltonian in the vicinity of the island whose 
eigenstates will approximate well those Floquet states which are localized in the vicinity of 
the island (see below). Simultaneously, time-periodicity of Floquet eigenstates assures that 
the exact Floquet states will not disperse. 

The existence of such wave packet Floquet eigenstates has been proven by an exact 
numerical diagonalization of the problem both for linear polarization (LP) |7[] and circular 
polarization (CP) |],|Hj of the microwaves. To allow their detection one should consider 
the ways of populating effectively such states. For a CP, where the direct optical excita- 
tion from a weakly perturbed low lying state is impossible (since the wavepacket is built 
from predominantly circular atomic states unaccessible from low lying states due to dipole 
selection rules) one should first prepare the atomic circular state [p2| , p3[ and then switch 
on the microwaves sufficiently fast ||. It has been shown that the wavepacket states may 
be populated in this way with about 90% efficiency P3fl . For a LP case, an addition of 
a static electric field allows to control the trajectory of the wavepacket in particular 



wavepackets moving along elongated, low angular momentum trajectories may be created. 
Such wavepackets may be accessible to a direct excitation from low lying states. 

The next important question is the possible mechanism of the detection of these states. 
At least three possible ways suggest themselves. Two of them utilize the residual decay of 
wavepackets, either via the spontaneous emission (treated both for LP |25| and CP [^6],^TJ) 
or the ionization. The former may not be efficient, since, at least for wavepackets moving 
on circular trajectories, the corresponding spontaneous emission rates are quite small ||2T] . 
In the ionization experiment, the population of the wavepacket state may be detected by 
a strong decrease in the ionization yield (since ionization rates of wavepacket states are 
typically very small [[||TU|). On the other hand, these rates fluctuate strongly (the mechanism 



of their ionization, via "chaos assisted tunneling" is discussed in detail elsewhere |19| , P0| ) - 
this may make their detection in the ionization yield quite ambiguous. 



By far the most promising method is the Floquet spectroscopy |27j i.e. probing, by a 
second weak microwave field, the structure of Floquet (dressed || by microwaves) states. To 
this end a precise estimation of the quasienergies of wavepacket states is necessary. An exact 
diagonalization of the problem gives all the Floquet states and a time consuming inspection 
of individual eigenvectors is necessary to identify the wave packet states. This process may 
be optimized by calculating properties of matrix elements of appropriately chosen operators 
but certainly it is desirable to have good semiclassical predictions for the quasienergies. 
For CP case those are given, to a very good accuracy from the harmonic approximation 
Hamiltonian [§-T0], this approach being, however, restricted to this particular system. 



For a general case of periodically driven system there is no simple unitary transformation 
which removes the time-dependence (as it is in the CP case) and the correct approach is to use 
approximate resonant Hamiltonians. The semiclassical quantization of such a Hamiltonian 
gives not only the good estimate for wavepacket states but allows for the classification of 
resonant states for systems of more than one dimensions. Recently using such an approach 
we could discuss the resonant dynamics in a realistic three dimensional (3D) model of a 
hydrogen atom in the LP case ||28|| . Similarly, we have discussed the control of wavepacket 



trajectories using an additional static electric field 

Up till now the discussion of nonspreading wavepacket states has been restricted to 
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linear and circular polarization cases only. The aim of this paper is to treat a resonant 
dynamics of a hydrogen atom, both semiclassically and quantum-mechanically in a general 
case of an elliptical polarization (EP). Apart from generalizing the notion of nonspreading 
wavepackets to an arbitrary EP, we discuss the full dynamics of quasienergies as a function 
of the ellipticity of the microwaves for a resonant case, being stimulated by recent ionization 
experiments ||29|| . Unfortunately, the experiments do not allow for a full selection of the initial 
state of an atom (states with all possible angular momenta are simultaneously excited) but 
this situation may improve in the near future. 

The EP case is highly nontrivial. For LP microwaves the conservation of the angular 
momentum projection onto the polarization axis, L z , makes the dynamics effectively two- 
dimensional (2D). For the CP case while L z is not conserved, the transformation to the frame 
rotating with the microwave frequency removes the explicit oscillatory time-dependence (see 
e.g. [j30| -|33|| ) ■ Both these simplifications are no longer possible in the general EP microwave 
field. In effect, the exact quantum diagonalization approach for the part of the spectrum 
corresponding to the strongly perturbed atomic Rydberg spectrum would require very big 
computer memory. Let us mention also that in the effective 2D LP case, the state of the art 
computations consider initial atomic states with the principal quantum number of the 
order of 20. 

For that reason we shall consider not the realistic fully 3D model of an atom but rather 
the restricted 2D model in which the electronic motion is restricted to the polarization plane. 
Study of such simplified models have been most successful in the past both for LP (where 
one dimensional model has been a main source of quantum results for a long time [p^,^5|,|6|) 
and in CP [|30|-|3"3"| where also the 2D, polarization plane restricted model has been utilized. 
Additional argument favoring the 2D model comes from our classical study of dynamics in 
EP microwaves |3B],|37| where comparison of 2D [^] and 3D analysis shows the similarity 
of physical phenomena in both cases. Simply put the perturbation is most effective if the 
polarization plane coincides with the plane of Kepler motion. 



II. THE SEMICLASSICAL VERSUS QUANTUM APPROACHES 

The Hamiltonian of the hydrogen 2D model atom driven by an elliptically polarized 
electromagnetic field reads in the dipole approximation and in the length gauge (in atomic 
units) 

p 2 + p 2 1 

H=— h F(xcosujt + ay smut), (2.1) 

2 r 

where r = \/x 2 + y' 2 while F and u denote the amplitude and the frequency of the mi- 
crowave field, respectively, a defines the ellipticity of the microwaves with a = (a = 1) 
corresponding to a LP (CP) limiting case. 

Using the Floquet theorem the solution of the quantum problem is equivalent to diago- 
nalizing the Floquet Hamiltonian 



H-i—)ip n = H F ij n = E n ij n (2.2) 
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with E n being the quasienergies while ip n time-periodic Floquet eigenstates. 

The details of the numerical method are described in the Appendix. In short, the calcu- 
lations proceed by expressing the Floquet eigenvalue equation in the scaled semi-parabolic 
variables 



u 2 — v 2 



x = A — - — , y = Auv. (2.3) 

where A is an arbitrary scaling factor. This allows to remove the Coulomb singularity and 
cast the Schrodinger equation into the generalized eigenvalue problem for polynomial like 
operators. Standard harmonic oscillator creation and annihilation operators allow for simple 
evaluation of matrix elements. The approach closely resembles that for the 2D atom in the 
CP field |33| except that the explicit time dependence is treated by the Fourier expansion. 



It is worth stressing that using the complex scaling parameter A one effectively realizes 
the complex rotation in the system which enables the exact treatment of the coupling to 
the continuum (ionization), for details see [0,^3J for LP and CP cases, respectively. In this 
paper, however, we shall concentrate on the Floquet level dynamics only for the sake of the 
comparison with the semiclassics. The analysis of the ionization phenomenon is left for a 
future publication. 

The semiclassical quantization of resonant dynamics closely resembles the similar proce- 
dure applied by us recently for the LP case fli8| , |28f . That in turn originates from a general 



prescription for EBK quantization of Floquet spectra p8| 

Starting with the Hamiltonian ( |2.1| ) we remove the explicit time dependence by going to 
the extended phase space ||. Defining the momentum p t conjugate to t (time) variable we 
get the new Hamiltonian 

n = H + p t (2.4) 

which is conserved during the motion. The quasi-energies of the system will be then the 
quantized values of 7i. 

As the next step we express the Hamiltonian in action-angle variables of the unperturbed 
Coulomb problem [|36|| . For the 2D model atom those are e.g. the canonically conjugate pairs 
(J, 9) and (L, <fi). J is the principal action (corresponding to the principal quantum number, 
n ). The corresponding angle, 9, determines the position of the electron on its elliptic 
trajectory. L is the angular momentum (equal to L z for the 2D motion in the x — y plane) 
while is the conjugate angle (the angle between the Runge-Lenz vector and the X~£LX1S ; 1.6. 
the main axis of the polarization ellipse). 

We shall consider below the case of the resonant driving, i.e. when the frequency of 
the Kepler motion u>x = 1/J 3 is close to the microwave driving frequency to. Applying 
the secular perturbation theory Q to average over the nonresonant terms one obtains an 
approximate resonant Hamiltonian (in the frame rotating together with the electron) of the 
form 

Hr = ~2J 2 ~ UJ+ Fr(L ' ^ a) COS ^ " ^ + Pt (2 ' 5) 

where 9 = 9 — ut while p t = Pt + ojJ . 7i T yields the pendulum like principal action motion 
with the strength and the equilibrium position determined by T(L, <p\ a) and 5 = 5(L, 0; a) 
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respectively. The pendulum Hamiltonian is obtained by additionally expanding 7i r around 
the center of the resonance island given by J = to^ 1 ^ 3 up to the second order in AJ = 
J — u; -1 / 3 but we do not apply this expansion. Both T and S depend on the initial shape 
and orientation of the electronic ellipse (via L, 0) as well as on the ellipticity of microwaves, 
a, and are given by p6| , |39 



r(L,<p;a) 



1 -±^V 1 ) 2 + 2cos2<p 1 -±^V 1 1 - 



a 



1 



21 V2 



(2.6) 



. (1 - a)V- X - (1 + a)V x 1 
(1 - a)V-i + (1 + a)Vi 



where 



VT^ 2 



V ±1 (J,L) = u- 2 ' 3 [ji{e) ± —^(e)]. (2.8) 
J7i(x) and j[{x) denote the ordinary Bessel function and its derivative, respectively, while 



e = yl — L 2 / J 2 = yl — L 2 u 2 / 3 stands for an eccentricity of an electronic ellipse. For 
completness let us mention that a similar approximation on the purely quantum level, leading 
to Mathieu equation has been performed for the CP case only in j|0] . 



The semiclasical quantization of (|2.5| ) is straightforward and follows closely the procedure 



described in detail elsewhere [28] for arbitrary m : 1 resonance. For 1 : 1 resonance considered 
here, the trivial quantization of p t , exploring the time periodicity of the system, yields 
additive terms ku to quasi-energies (different values of k correspond to different Floquet 
zones). As discussed first in see also p8| , the orbital motion in (J, 6) variables (along 
the perturbed Kepler ellipse) is much faster than the modification of the ellipse shape and 
its movement (precession) as described by the motion in (L,<ft). Thus the slow and fast 
motions may be adiabatically decoupled. Making Born-Oppenheimer approximation and 
using standard WKB rules, the fast (J, 9) motion is quantized taking the Maslov index 
v = 2 (corresponding to librations, i.e., we quantize states inside the resonance island). 
Being interested in resonantly localized states we shall consider later the ground state of the 
radial motion, only. For the slow angular (L, 0) motion we take the Maslov index \x = or 
2 for a rotational or librational motion, respectively. 

Similarly, as in the LP microwaves |p8[| , it is easier to quantize first the slow motion 
generated by constant values of T(L, (f>; a) and later treat the fast motion. Such a procedure 
is justified since quantizing the fast motion one takes T(L, <f); a) as a constant quantity and 
thus the order of the quantization does not matter. 

The existence of the resonance island in (J, 9) space ensures that the radial motion 
is localized. So the remaining analysis should concern the angular (L, <fi) motion which 
reflects the slow evolution of the Kepler ellipse. The structure of the (L, 0) space influences 
values of quasi-energies as well as the structure of corresponding semiclassical eigenstates. 
Trajectories in the (L, <p) space are determined by constant values of T also responsible for 
the size of the resonance island in (J, 9) space (recall that the island's size is determined 
approximately by y/FT). 
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Before presenting the results let us define scaled variables, typically used as a convenient 
parameterization of the microwave ionization problem (since the dynamics scales classically 
pT2|). The common choice is to link the scaling to the initial energy of the electron ||35|| . 



For an unperturbed 2D hydrogen atom the eigenenergies are given by E n = — l/2(n+ 1/2) 2 
with n being non-negative integer. This allows us to define scaled variables via the principal 
quantum number of the initial state hq via 

lu = uj(n + 1/2) 3 (2.9) 
F = F(n + 1/2) 4 (2.10) 
L = L/(n + 1/2). (2.11) 

In particular note that the scaled angular momentum Lq may take values from [—1,1] interval 
with extremal values corresponding to circular orbits on which the electron moves in two 
opposite directions. 



III. RESULTS AND DISCUSSION 

To compare the semiclassical predictions to the exact quantum calculations, we consider 
the n = 21 manifold of our 2D model atom. This value is a compromise between the 
requirement to be in the semiclassical, large n regime and computer memory limitations 
(size of the Floquet matrix to be diagonalized) . For resonant driving we take the microwave 
frequency to be u = ujk = l/(^o + 1/2) 3 ) (^o = !)• 

In Fig. m we present values of V as a function of the scaled angular momentum, L 
and the <fi angle for two different values of the degree of the field ellipticity given by a. 
Equipotential curves of V corresponding to semiclassical states originating from n = 21 
hydrogenic manifold are also shown in the Figure. Due to the specific form of the resonance 
Hamiltonian ( |2.5| ) those lines are independent of the microwave amplitude, F. 

For a = the (L, <ft) space is symmetric with respect to the L = axis since, in the LP 
case, dynamics is not affected by the direction of the rotation of an electron. The left column 
in Fig. m presents the results for a = 0.1, i.e., the case very close to the LP problem. Note 
the presence of four stable fixed points: two of them, L ~ and = ir/2, 3ti/2, correspond 
to almost straight line orbits oriented perpendicularly to the main axis of the polarization 
ellipse, the other two fixed points correspond to circular orbits, \Lq\ ~ 1 (for such orbits 
is a dummy variable), with an electron rotating in the same or opposite direction to the 
direction of the rotation of the field vector. In the vicinity of these fixed points the electronic 
motion is "shape" localized (since the motion remains in the vicinity of a given (Lq, 0) point 

the eccentricity of the orbit, e = \Jl — Lq, as well as its orientation with respect to the 
polarization ellipse, given by is approximately preserved). Whether the motion is well 
localized in its radial motion in (J, 9) variables, i.e. along the ellipse depends on the size 
of the resonance island, given by V FT as mentioned above. If the island in (J, 9) space is 
sufficiently large (comparable in size to h or larger) the quantum state will be localized in 
the island. On the contrary, too small resonance island cannot lead to localization in (J, 9) 
space, possible quantum states will spread over the whole ellipse (then also the semiclassical 
approach used becomes obviously non adequate). Thus to observe interesting, nonspreading 
wavepackets the resonance island must be large enough. 
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In this respect, two fixed points corresponding for small a to almost circular orbits, |Lo| ~ 
1 lay at the maxima of T (compare the left top panel in the figure) and for sufficiently large 
F may enable a strong radial localization. The other two stable fixed points, corresponding 
to elongated orbits lay at local minima of T. There exist also two unstable fixed points 
around Lq « and = 0, tt, they form the origin of the separatrix dividing the space. 

When a increases the fixed points situated initially around Lq « and <fi = ir/2, 37r/2 
move in the direction of greater negative values of Lq (see the right column in Fig. p, thus, 
the eccentricity of the corresponding orbits decreases. The circular orbits do not change the 
shape but the resonance island in (J, 9) space associated with L ~ 1 becomes larger while 
that associated with L ~ — 1 becomes smaller as reflected by the greater and smaller values 
of T, respectively. Note also that the islands in the (L, 0) space containing librational states 
shrink with the increase of a, hence, they can support less and less semiclassical states - 
while a increases librational states vault over the separatrices and become rotational. In the 
limit of a = 1 the islands disappear and there exist only rotational states. 

Consider the level dynamics corresponding to the change of a. Fig. |2| shows the level 
dynamics for the group of states originating from the tiq = 21 hydrogenic manifold. We 
take the scaled field amplitude Fq = 0.03, which for the CP case corresponds already to a 



significant ionization yield |3S| and may be, therefore, considered as a quite large value. For 
the LP case, i.e. a = 0, all states are doubly degenerate. Those of them which correspond 
to rotational states in the (L, <p) space are degenerate because the change of the rotation 
of an electron does not affect dynamics in the LP problem [the (L, <fi) space is symmetric 
with respect to the Lq = axis]. The remaining librational states are degenerate as the 
two islands, around L = 0, = 7r/2 and Lq = 0, <fi = Stt/2, are identical and support 
identical states. Of course the degeneracies may be removed due to tunneling processes, 
for instance, the wavepackets localized on circular orbits (the highest levels) rotate in the 
opposite directions and belong to distinct semiclassical states but quantum mechanical states 
are the symmetric or antisymmetric linear combination of them. Thus two wavepackets 
propagating on a circular orbit in the opposite directions correspond to a single eigenstate 
(these wavepackets were already discussed in [43] for the limiting LP case only). Similarly 
tunneling effects affect the librational states, i.e. quantum states are the linear combinations 
of solutions in each island. 

An increase of a removes the semiclassical degeneracy of rotational states as one can see 
in Fig. |2| - observe splittings of upper levels in the manifold. It corresponds to the broken 
symmetry with respect to change of a direction of rotation of an electron. Quasi-energies of 
states corresponding to the motion in the opposite sense to the rotation of the field vectors, 
i.e. corresponding to negative values of L , move down while those with positive L go up. 
It is a consequence of the behavior of T, i.e. the greater T the greater the quasi-energy value, 
see Figs. [I]. Observe also that the degeneracy of all librational states (lower levels in Fig. ^) 
is not immediately removed after a slight change of the field ellipticity from a = but is 
removed successively during an increase of a. This reflects the shrinking of the corresponding 
islands (the situation mentioned above) which causes the librational states to vault over the 
separatrix and become rotational. The levels with the smallest energy difference correspond 
to the librational and rotational orbits closest to the separatrix. The narrowing of the level 
spacing in their vicinity is just a consequence of the slowing down of the classical motion 
In the limit of a = 1 there is no degeneracy in the manifold. 
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Note that the energy splitting of the manifold is the largest in CP case while the smallest 
for LP. It is associated with the corresponding strength of the perturbation and simply 
expresses the dependence of T on a. 

Exact quantum results coming from the numerical diagonalization of the Floquet Hamil- 
tonian are presented also in Fig. ^. Among the multitude of Floquet states appearing in the 
same energy range only those with the largest overlap on the initial manifold are plotted 
for clarity. One can see that the agreement with the semiclassical predictions is very good 
except in the region of broad avoided crossings (with other levels - partners in the crossing 
eliminated by the overlap selection) which appear in the upper part of the Figure. The 
semiclassical method does not take into account the interaction of the considered n = 21 
manifold with other manifolds (the method describes only a single resonant manifold), thus 
such avoided crossings have no chance to appear in our semiclassical calculations. For a 
more accurate comparison we plotted, in Fig. [| quasi-energies for two different values of 
a separately. The agreement between the semiclassical and quantum results is very good 
for high lying levels (greater values of T) while for lower levels some differences appear as 
expected due to small size of the resonance island. 

In the limiting case a = 1, for CP case, another analysis of resonant dynamics is pos- 
sible, namely a harmonic approximation around a stable equilibrium point in the rotating 
frame || (corresponding to the expansion about a stable fixed point Lq — 1, <fi arbitrary 
in our picture). It is interesting to compare both approaches. The present semiclassical 
quantization of the resonance hamiltonian is certainly a more general approach valid for a 
general EP case and for the whole manifold. The harmonic expansion is limited to the CP 
case and valid for highest lying states in the manifold only. On the other hand this expan- 
sion is quadratic in deviations from the equilibrium point but nonlinear in the microwave 
field F strength while the resonance Hamiltonian approach is a first order in F expansion. 
Clearly in the deep semiclassical regime (large no) and for sufficiently large F$ the harmonic 
expansion approach yields a better approximation for wavepacket states in the CP than the 
present resonant Hamiltonian analysis (yielding, however, little information on the energies 
of other states in the resonant manifold). On the other hand, quite surprisingly, we have 
found, by a direct comparison of numerical values, that for F around 0.01 or 0.06 and no up 
to 30-40 the resonant Hamiltonian yields semiclassical values closer to exact quantum data 
than the harmonic expansion. Thus for intermediate range of no and Fq values the resonant 
Hamiltonian quantization is surprisingly good in accuracy, yielding, at the same time, the 
predictions for the whole resonant manifold and for an arbitrary polarization. 

To complete the picture we show in Fig. ^ the level dynamics of the same no = 21 
manifold versus scaled field amplitude, F , for a = 0.6. Quantum results are presented 
together with the semiclassical ones. Again, for the resonance island in orbital (J, 9) motion 
sufficiently large to capture quantum states, semiclassical results reproduce quantum ones 
quite well even beyond a classical chaos border (of course except avoided crossings). 

Exact quantum Floquet matrix diagonalizations, performed using the Lanczos code yields 
not only the eigenvalues but also the corresponding eigenvectors. Their time evolution may 
be visualized j|5| to confirm directly that indeed the Floquet states localized both in (J, 9) 
and in (L, <ft) spaces correspond to the nonspreading localized wavepackets. Consider again 
the level dynamics presented in Fig. |2|. As discussed above the higher lying state in the 
manifold corresponds to the wavepacket motion of the electron on the circular orbit in the 
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direction coinciding with the direction of rotation of the electric field vector. And indeed 
such a motion is revealed by the plots of the corresponding Floquet state (see Fig. [|, upper 
row). This Floquet state changes little with change of the ellipticity of the microwaves, a, 
and in the limiting case of CP becomes a well known CP nonspreading wavepacket |||9| . In 
the opposite limit of a = i.e. the LP microwave case this state becomes almost degenerate 
with another one corresponding to the different direction of the electron rotation. The two, 
almost degenerate exact Floquet states are linear combinations of the two wavepackets (at 
least in 2D) as discussed in . 



It seems interesting to see what happens to the second member of the pair as polarization 
is changed from being linear to elliptical (a increases). The two states separate fast in a 
and each of them represents a distinct motion. As mentioned before the state going down 
in energy corresponds to wavepacket moving on the circular orbit in the direction opposite 
to the field. It undergoes a series of avoided crossings with other states of the manifold 
loosing progressively its localized character. Still, for a not too large, and far from avoided 
crossings its wavepacket character is clearly visible (compare lower row in Fig. [5]). This 
state looses its wavepacket character when the librational islands (compare Fig. [I]) move 
sufficiently far down so as rotational states with large negative Lq disappear since there is 
no "space left" for them due to the finite value of h. As may be seen by comparison of 
Fig. [I] and Fig. [|this transition occurs around a = 0.45 for n Q = 21, manifesting itself in the 
quantum data (Fig. §) by relatively broad avoided crossings encountered by the wavepacket 
state. The transition point is, of course, h dependent, for smaller h (larger n ) the rotational, 
wavepacket-like, state with large negative Lq exists for larger a. It is worth stressing that for 
a below this critical, no (and thus driving frequency) dependent value there are two distinct 
Floquet states corresponding to two wavepackets moving in opposite directions. 



IV. CONCLUSIONS 

To conclude we have shown, by the comparison with the quantum numerical results, 
that the proposed semiclassical method, based on the first order resonant Hamiltonian gives 
good quantitative predictions. Numerical calculations produce nothing but numbers and 
an advantage of a semiclassical analysis is an understanding what physics is hidden behind 
them. Apparently the behavior of H atom placed in EP microwave field, in the range of 
parameters used above, is determined by the underlying classical dynamics, especially, one 
may build wavepackets which follow classical trajectories without spreading. We have shown 
that the wavepacket discovered in CP microwaves |||9| corresponding to the motion of the 
electron on a circular orbit in the same direction as the rotation of the microwave field 
exists also in EP case. For a sufficiently small ellipticity of the microwave polarization, a, 
there exist also another wavepacket state, corresponding to electron rotating in the opposite 
direction to the microwave field. In the limiting LP case both wavepackets coalesce, the 
exact Floquet states correspond semiclassically to linear combinations of two wavepackets 



propagating in the opposite direction on the circular orbit 

The analysis presented is restricted to the two-dimensional model atom (similar assump- 
tion is implicit in P5|), its validity for the real three-dimensional atom is an open question. 
Certainly, in the limiting LP case, due to rotational symmetry with respect to the field 
axis, the two wavepackets moving in the opposite directions in 2D lead in 3D to a doughnut 
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shaped localized function oscillating between north and south pole of the sphere (assuming 
a vertical polarization of LP microwaves) |2£| . For CP, on the other hand, the wavepacket 



motion remains essentially two-dimensional ||10||. It will be most interesting to see how the 



third dimension affects the dynamics also for a general case of EP microwave field since 
already classical studies p6| , |37| indicate that some qualitative differences may appear. This 
subject is left for future studies. 

Still the present 2D analysis shows that the nonspreading wavepacket states localized on 
circular orbits are not restricted to circular polarization of microwaves only. Thus the perfect 
circular polarization is not essential for a possible experimental observation of nonspreading 
wavepackets. 

Interestingly, in the EP case, there appears a possibility of the angular localization in 
minima of the effective potential T with the position of the minima being a dependent 
(compare Fig. [I]). If the creation of localized wavepackets in such minima were possible it 
would allow to control the shape of the trajectories on which the wavepackets propagate 
by a change of microwave polarization, and not by an additional static field as proposed in 
||18|| . Study of the corresponding wavefunctions indicates, however, that due to the small 
resonance island width the localization along the ellipse [i.e. in (J, 8) space] is not very 
effective at least in the range of n and F studied by us. 
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VI. APPENDIX 

Consider the Floquet Hamiltonian, Hp, as defined in ( f2.2| ) with the hamiltonian given by 
( 12.1|) . For the sake of efficiency of the numerical diagonalization it is convenient to rewrite the 
Floquet Hamiltonian in the velocity gauge by applying a standard unitary transformation 

urn 



d pl+pl 1 F r F 2 , 

H F = -i—-\ — y - 1 \p x $mu)t -p y a cos ut \ + (a +1). (6.1) 

at 2 r uj Auj 1 



Next we pass to scaled semiparabolic coordinates ( |2.3[ ) which allows us to express the Floquet 
eigenvalue equation 

H F | u n (t)) = E n | u n {t)) (6.2) 
clS db generalized eigenvalue problem 

(A-EJ3) | u n {t)) = (6.3) 
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where 



A= (u 2 + v 2 )H F and B = u 2 + v 2 



(6.4) 



with 



pI + pI 



2 . d 



2A 2 (u 2 + v 2 ) A(u 2 + v 2 ) dt Alu(u 2 + v 2 



■[(up u — vp v ) smut — a{vp u + up v ) cosut] 

(6.5) 



(with the additive pondoromotive term omitted). 

Floquet states are periodic in time (with the period T = 2tt/u) thus we may expand 
| u n {t)) in the Fourier series 



-iKut 



(6.6) 



This allows to cast Eq. ( |6.3|) into the equivalent set of coupled equations 



Pj+Pj 2 
2A 2 A 



where 



- Kuj(u 2 + v 2 



<) + (W s - W c ) | - (W. + W c ) I = £^(« a + v 



(6.7) 



-{up u - Vp v ) 



2coA 

F . 

a^j^{vp u + up v ) 



(6.8) 



All terms in the above equation have the polynomial form in coordinates and momenta. This 
suggests to use the harmonic oscillator basis for an efficient evaluation of matrix elements 
39| . The method becomes then analogous to the treatment of the circular polarization 
case discussed in detail elsewhere |33|. The resulting generalized eigenvalue equation is 



diagonalized using the Lanczos code which allows for an extraction of eigenvalues in a selected 
energy range and the corresponding eigenvectors. For completeness let us mention only that 
the Floquet Hamiltonian is invariant under the generalized parity transformation, i.e. the 
parity transformation combined with the translation in time by ir/u. Thus A and B matrices 
may be split into uncoupled matrices that are two times smaller. This makes the numerical 
calculations more efficient. 
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FIGURES 



FIG. 1. Two dimensional hydrogen atom illuminated by resonant, uq = 1, elliptically polarized 
microwaves. The effective scaled perturbation Tq = V /{jiq + 1/2) 2 is plotted as a function of the 
scaled angular momentum, Lq, and the angle, <ft, between the Runge-Lenz vector and the main 
axis of the polarization ellipse - upper row. Bottom row shows equipotential curves of the angular 
part r of the Hamiltonian 7i r , Eq. (|2.5| ), representing the slow evolution of the Kepler ellipse. 
The curves correspond to semiclassical states originating from the tiq = 21 hydrogenic manifold. 
Columns correspond to the different ellipticity of the microwaves a = 0.1 (left) and a = 0.6 (right). 

FIG. 2. Two dimensional hydrogen atom driven by resonant, loq = 1, elliptically polarized 
microwaves. Level dynamics, versus a (i.e. the degree of the field ellipticity), of the semiclassi- 
cal quasi-energies (full lines) of the states originating from the tiq = 21 hydrogenic manifold for 
Fq = 0.03 compared with the exact quantum results (dotted lines). 

FIG. 3. Comparison of the semiclassical quasi-energies (circles) originating from the unper- 
turbed uq = 21 manifold for Fq = 0.03 and ujq = 1 with the exact quantum values (crosses), at 
different values of the degree of the field ellipticity a = 0.1 (a), 0.9 (b). Integer index p counts 
consecutive states in the perturbed manifold. 

FIG. 4. Level dynamics of the exact quantum quasi-energies (dotted lines) in the vicinity of 
the resonant manifold emerging from no = 21, compared with the semiclassical prediction (full 
lines), for Fq = . . . 0.06 and a = 0.6, loq = 1. Note that the maximum field amplitude clearly 
exceeds typical ionization thresholds measured in current experiments at the principal resonance. 
The semiclassical prediction accurately tracks the exact solution across a large number of avoided 
crossings. 

FIG. 5. Nonspreading wavepackets for 2D hydrogen atom illuminated by the elliptically po- 
larized microwaves with the amplitude Fq = 0.03, frequency luq = 1 for uq = 21 and ellipticity 
a = 0.4. Top row - the exact Floquet state corresponding to a nonspreading wavepacket moving 
on a circular orbit in the direction of the rotation of the microwave field at times cot = 0, 7r/4, tt/2 
from left to right. This wavepacket corresponds to the one known for the circular polarization and 
may be obtained from the latter by a change of the microwave polarization. Bottom row - another 
exact Floquet state corresponding to a nonspreading wavepacket moving in the opposite direction 
(shown at the same times ut = 0,7r/4, it/2 from left to right). Note that while the former almost 
preserves its shape in temporal evolution, the latter becomes significantly distorted. Still being 
an exact Floquet time-periodic eigenstate it regains its shape (as depicted e.g. in the bottom left 
corner of the figure) every period of the microwave. The size of each box is ±800 Bohr radii in 
both x and y directions. 
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